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Abstract 


Pathways-reduced analysis is one of the techniques used by the Fispact-II nuclear 
activation and transmutation software to study the sensitivity of the computed inven¬ 
tories to uncertainty in reaction cross-sections. Although deciding which pathways are 
most important is very helpful in for example determining which nuclear data would 
benefit from further refinement, pathways-reduced analysis need not necessarily de¬ 
fine the most critical reaction, since one reaction may contribute to several different 
pathways. This work examines three different techniques for ranking reactions in their 
order of importance for determining the final inventory, comparing the pathways-based 
metric (PBM), the direct method and one based on the Pearson correlation coefficient. 
Reasons why the PBM is to be preferred are presented. 
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1 Introduction 

Fispact-II is a software suite for the analysis of nuclear activation and transmutation events 
of all kinds [TJ. The present work focuses on its use in sensitivity studies of nuclear inventory 
calculations; these employ the Bateman model for the evolution of the inventory of a target 
subject to irradiation by an imposed flux of projectile particles, always neutrons herein. In 
ref [2] it was established that the pathways-reduced approach jam to such studies, almost 
invariably gives very close agreement with Monte-Carlo sensitivities computed using full 
Bateman, i.e. accounting for all nuclides and pathways. Pathways-reduced models are, 
following Eastwood and Morgan [3J, identified by a graph-based approach which determines 
the key reaction pathways determining the inventory at a given time and eliminates from 
consideration those nuclides which do not lie on this reduced set of pathways. 

The pathways-reduced metric is a sensitivity method in the respect that implicitly it se¬ 
lects a set of the most important nuclide reactions. A wide range of other sensitivity methods 
has been used by the nuclear industry as shown by literature reviews by Helton et al [5], 
see also Cacuci and Ionescu-Bujor PIT). Indeed, sensitivity analyses are available as part 
of nuclear industry software packages such as for example DAKOTA j8j and SCALE |9j. 
General purpose sensitivity software is also available, such as OpenCossan [ 10] . 
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A key input to most techniques considered herein is an estimate of the uncertainty in the 
reaction cross-section. The determination of such uncertainties is a challenging subject in its 
own right, hence it is important to examine techniques that can identify which reactions most 
require further examination. The present work represents a comparison of three different 
techniques that exploit the pathways-based reduction for the nuclear activation problem. 

Fispact-II can access uncertainty data, typically the standard deviation, for the vast 
majority of reaction cross-sections in the EASY-II database m however no information is 
currently passed concerning pure decay reactions. This reflects the fact that half-lives are 
often very accurately known. There are other reactions in the database for which a value of 
zero uncertainty is found, usually indicating that no information is available. There are thus 


difficulties in making the comparison, the implications of which are discussed in Section 2.4 


To proceed further with this introduction, it is efficient to introduce the time evolution 
(rate or Bateman equation) for a nuclear inventory X 


dX 

dt 


AX 


( 1 ) 


where X is the vector of nuclide numbers, and A is the matrix of nuclear interaction coef¬ 
ficients for both induced reactions and spontaneous decays. Note that one coefficient A t j 
of A may represent several different nuclear reactions, since the equation involves an average 
over a spectrum of energies (of neutrons in the present work, although other elementary 
particles may be considered in general). Hence the term ‘interaction’ is used to cover all 
effects generating nuclide i as the child of parent j. It is worth noting that although i pre¬ 
cedes j alphabetically, reactions throughout this work will except for the A^ be described 
in parent-child order. In general, the coefficients Aij may change with time as the incident 
neutron flux changes. 

All the techniques for ranking the interactions A^ are most easily understood in the con¬ 
text of a single constant irradiation in the time interval (0 ,tf), producing an inventory X(tf). 
Different aspects of the inventory, such as heat production or kerma, may be studied using 
FlSPACT-II, but for illustrative purposes it is sufficient to consider only the total activity 


Q — ^ A fc AA(t/) 

k 


( 2 ) 


where A k is the decay rate of the nuclide X*,; A k is zero for stable nuclides and A*, = log e 2/r^ 
for unstable ones, where is the half-life. Although attention focuses here on FlSPACT-II, 
the Bateman approach is standard in that most packages with a claim to generality, not 
only DAKOTA and SCALE in the U.S. but also ANSWERS [12] with FISPIN in the U.K., 
include solvers for the problem. The pathways-based analysis technique studied here could 
in principle be implemented in any of these codes. 

The three different ranking techniques are described in the next Section [2} Apart from the 
use of pathways-based reduction, there is novelty in the calculation of the direct sensitivity, in 
that the matrix Frechet derivative is used in its computation, see Appendix, rather than the 
more usual decoupled direct method DDM of Dunker im The application of the techniques 
to the wide range of test cases first introduced in ref [2] is illustrated in Section [3j Lastly 
Section [4] compares the utility of the different techniques. 
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2 Sensitivity Measures 

2.1 Pathways Based Metric 

The Pathways Based Metric (PBM) is calculated quite simply from the output of the 
pathways-reduced approach, which includes a listing of each pathway and its percentage 
contribution to the active nuclide at its termination. For a given interaction A t j , all the 
number N p of pathways upon which it lies are identified and the PBM calculated as 

N p 

SpBM = E PiXtX t hi (3) 

k= 1 

where pi is the fractional contribution of pathway l to the number of atoms X t (evaluated at 
time tf ) in the inventory with decay rate A t and the indicator matrix I^i = 1 or 0 depending 
whether or not a reaction contributing to the interaction lies on the pathway. 

Figure [l] illustrates how the definition works in a simplified case where irradiation of an 
initial sample consisting of X± atoms of nuclide 1 and X 2 atoms of nuclide 2 produces an 
inventory containing numbers X 5 and Xq of radioactive nuclides 5 and 6 respectively, with 

3 important pathways. The first pathway contributes piX 3 atoms and the third p 3 X 3 atoms 
of nuclide 5. (Supposing that other pathways are unimportant, pi + p 3 ~ 1 and p 2 ~ 1.) 
The sensitivity of the inventory to the reaction with coefficient A 32 (large arrowheads in 
Figure [Tj) , is then 

$ 1 ?bm = P1X5X5 + p 2 XeX 6 ( 4 ) 

where A 5 is the decay rate of nuclide 5 etc. 

This technique required special modification to FlSPACT-II for its implementation, which 
was facilitated by the object-oriented design of the Fortran-95 code. For the purposes of 
initial investigation, the loops which are identified by the graph-based approach used by 
FlSPACT-II are ignored. 

2.2 Direct Method 

The Direct Method (DM) works directly with the tensor describing the rate of variation of 
the nuclide X & with respect to nuclear reaction coefficients. For initial investigative purposes 
it is sufficient to consider the partial derivative with respect to A iy Differentiating Eq. ([Tj) 
with (z, j) regarded as fixed, gives 

d ( dX \ _ dX dA x 

dt \dAij) ~ A 3A ij + dAj ^ } 

If the sensitivity of the total activity is required, then using Eq. ([2]) , this is 

Am = E^ 9 X£ j1 (6) 

k u ^*j 
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Figure 1: Example illustrating how the Pathways Based Metric Spbm is calculated when 
the three pathways shown are highlighted by the pathways-reduced analysis. The pathways 
are numbered in order 1, 2 and 3 from the top. Pathway 1 starts with nuclide 2 and generates 
via a sequence of reactions involving nuclides 3 and 4, p\X 3 atoms of nuclide 5, whereas the 
other pathways begin with nuclide 1 and generate nuclides 6 and 5 respectively. The larger 
arrowheads indicate reactions with the coefficient of interest A 32 . 

In the decoupled direct method, Eq. (J5]) is solved for dXk/dAij using a method which exploits 
the sparseness of dA/dA^ = Sij in the present context. However, it is also possible to express 
Sbm bi terms of the matrix Frechet derivative as explained in the Appendix, viz. 

S$ DM (t f ) = t, £ Ey)X(O) (7) 

k 

where L exp is the matrix Frechet derivative as defined in the Appendix where E, ? is also 
defined. Eq. ([?]) defines the Frechet direct method. Similarly to the PBM, this technique 
required modification of FlSPACT-II to output the matrix A in a format suitable for input 
to MATLAB P3]. 

2.3 Pearson Derived Method 

The Pearson technique for ranking sensitivities starts with the definition of the Pearson 
product-moment correlation coefficient for a set of N s samples {(A S ,Q S ) : s — 1,..., N s }, 
viz. 

r = T,,{Qs-Q){A.-A) 

(N 8 - 1)AQAA { 1 
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where the suffix ij on r and A is to be understood, overbar denotes average and A denotes 
the standard deviation of the distribution so that for example 


A 



( 10 ) 


( 9 ) 


The coefficient r l3 is by definition always less than or equal to one, and a magnitude of r 
close to one indicates strong linear correlation. 

However, it is the proportionality constant corresponding to dQ/dAij that is of initial 
interest. It is worth cautioning that although the definition implicitly implies a linear rela¬ 
tion, there is no guarantee of this. However, in order to proceed, linearity will be assumed, 
viz. 


Q - Q = r(A - A) 


( 11 ) 


and substituting in Eq. (J8|, it follows that 



( 12 ) 


It follows that the output of the Monte-Carlo sensitivity calculations may be used to rank 
the different interactions by computing ry,-/AH (note that A Q is the same for all the A l3 vari¬ 
ations in the standard approach described in ref [2j). 

The calculation of the Pearson coefficient r is well-known to be sensitive to round-off 
error. To avoid modifying the software, the coefficient is computed using output values from 
FlSPACT-II given only to 6 significant figures by default. This accuracy is the maximum that 
can be expected from the numerical integration of the rate equation which is is constrained 
to an accuracy of one part in a million. It was found that splitting the separate contributions 
of A s and A to Eq. ([8]) led to unacceptable cancellation due to round-off effects (although it 
was verified that round-off was not a similar issue for Q s and Q ). 

2.4 Comments upon the Different Metrics 

The main distinction between the PBM and the other two measures is that the pathways- 
based method is ‘global’, capturing the whole variation of the inventory as parameters are 
varied, although having the disadvantage that it cannot measure sensitivity to diagonal 
entries of A. The other two techniques are more local, indeed the DM returns directly only 
a coefficient at the mean of the distribution of Q. The Pearson method is somewhere in- 
between, using global variations, but making a local linear assumption about the mean. This 
complicates the comparison in the next Section [3j 

The principal comment to be made concerning the comparison is that, corresponding to 
the lack of sensitivity to element Aij when it is zero due an absence of interaction between 
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nuclides i and j, a large sensitivity in the local sense may be inconsequential for the total 
activity Q if the corresponding is relatively very small. However, the two more local 
estimates (Eq. ([6]) and Eq. (12)) for dQ/dA^ should be directly comparable. 

Main interest attaches to global measures such as Spbm■ The FDM approach may be 
used to produce an equivalent ranking by scaling by the estimated error in the coefficient, 


viz. 


qij _ qij ( ( AL\ A 

FDS ^ TT'nA/T l _ I ' AA 


J FDM 


V100/ 


(13) 


where e %3 is the percentage error in the distribution of the coefficient A t j. FlSPACT-II returns 
and Aij by combining the uncertainties in the reaction coefficients corresponding 


both 
to A^. 

From Eq. (12), a ranking based on the Pearson coefficient r should also be comparable 
to Spbm, if it is scaled similarly, viz. 


S PR3 ~ Tij ■ (^q) 



(14) 


In practice it is found that Sp RS ~ ry,-. 

Note that for interactions for which no uncertainty information is available, a Pearson 
coefficient cannot be computed, nor is Spps useful. The coefficient Spbm ma y be non¬ 
zero, but this relies on the interaction’s lying on a pathway important for other reasons. 
Interactions without accompanying uncertainty information will therefore largely have to be 
ignored in this work. 


3 Sensitivity Calculations 

3.1 Details of Cases 

The test cases are taken from ref [p] and involve several different nuclide mixtures designed 
to be indicative of a wide range of activation problems, see Table [lj As indicated, all but 
one of the mixtures consisted of 1kg of material subject to a neutron flux of 10 15 cm _2 s _1 , 
for a year, without any cooling period. 

The mixtures are used in six test cases, with the Alloy case extended to include a cooling 
phase. Each test case is run using the full TENDL 2013 data from the EASY-11 database HU 
with pathways analysis to identify the important reactions, the numbers of which are listed 
in Table [2j As in ref [2j, Monte-Carlo solution of the reduced problem, investigating the 
distributions of the important reaction rates specified in the newer database, was then per¬ 
formed in the sequence of increasing sample size per reaction, N x = 10, 40, 160, ... up to 
the maximum value specified in the table. Indications from ref [2] and work which may be 
published elsewhere indicate that the pathways-reduced results agree to at least two (and 
often three) significant figures with those obtained by sampling the full problem, at less than 
a thousandth of the computational cost. As might be expected from the large maximum 
number of samples N s employed, the distributions of reaction rates actually sampled usually 
agree in the mean to 4 significant figures with the nominal database values. 
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Table 1: Test cases. Each consists of numbers of atoms of the listed elements with their 
natural abundances of nuclides, given as percentages by mass of the whole. 


Test 

Label 

Constituents of Mixture 

Sample 

Mass 

Irradiation 

Period 

Cooling 

Period 

Neutron flux 

—2 —1 
cm s 

Alloy 

Fe 40.0 : Ni 20.0 : Cr 20.0 : Mn 20.0 

1 kg 

lyr 

0 

10 15 

Alloy+c 

Fe 40.0 : Ni 20.0 : Cr 20.0 : Mn 20.0 

1 kg 

lyr 

lyr 

10 15 

Fe 

Fe 

1 kg 

2.5 yr 

0 

10 15 

LiMix 

Li 40.0 : Be 30.0 : 0 30.0 

1 kg 

lyr 

0 

10 15 

WMix 

W 20.0 : Re 20.0 : Ir 20.0 : Bi 20.0 
: Pb 20.0 

1 kg 

lyr 

0 

10 15 

Y203 

Y 78.74 : 0 21.26 

lg 

300 s 

0 

1.116 x 10 10 


Table 2: Test cases statistics. Monte-Carlo sampling by Fispact-II has a sample size 
determined by the number of reactions examined. 


Test 

Label 

/, Reactions 
Examined 

Matrix 
A Size 

Max. N x , Samples 
per Reaction 

N s , Total 
Sample 

Alloy 

84 

51 

640 

53 760 

Alloy+c 

50 

38 

640 

32 000 

Fe 

27 

24 

640 

17280 

LiMix 

17 

21 

640 

10 880 

WMix 

71 

63 

640 

45440 

Y203 

13 

16 

2 560 

33 280 





3.2 Results 


This section presents results for each of the test cases in turn, in the alphabetic order specified 
in Table [l] Attention is drawn to the fact that the Y203 case is the simplest in terms of 
pathways, and contains extra explanation. 

For each test case there is a table of sensitivity rankings ordered by Frechet derivative 
amplitude and a graph of rankings ordered by Spbm■ The table enables a larger range of 
interactions to be compared, since the graphs become hard to interpret once the number 
of plotted interactions exceeds about ten. Note the convention (except for the Y203 case) 
that all three methods must provide a ranking for the comparison to be plotted. So in the 
figures the ten highest-ranked cases plotted may include reactions significantly smaller in 
importance than the tenth. 

A general feature of all graphs comparing rankings by the different techniques is the 
symmetry about the mid-line labelled PBM. The appearance of “V” and “A” patterns in¬ 
dicates that although the more local measures may not agree with Spbm, they do themselves 
correlate well. 

For three of the test cases, Alloy+c in Section 3.2. 2[ WMix in Section 3.2.5 and Y203 
in Section 3.2.6 further results of analysis are presented to help understand the effect of 
sampling and round-off on the calculation of Spps■ I n addition, a table of sensitivity rankings 
ordered by Spbm and a graph of rankings ordered by Frechet derivative also appear in 


these two sections. (This information is omitted from the other four sections Section 3.2.1 
Section 3.2.3, Section 3.2.4 and Section 3.2.6 to save space.) 


3.2.1 Alloy 

See Table [3] and Figure [2j 
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Table 3: Alloy case. Rankings for different methods. 




Sensitivity 

Parent 

Child 

FDS 

PBM 

PRS 

Fe-56 

Mn-56 

1 

5 

1 

Ni-58 

Ni-57 

21 

37 

57 

Mn-55 

Cr-55 

10 

11 

10 

Mn-55 

V-52 

12 

18 

12 

Mn-55 

Mn-56 

9 

1 

9 

Cr-52 

Cr-51 

17 

19 

24 

Cr-52 

V-52 

5 

10 

5 

Ni-58 

Co-58m 

3 

4 

3 

Mn-55 

Mn-54 

11 

6 

11 

Fe-56 

Fe-55 

14 

15 

15 

Ni-58 

Co-57 

8 

9 

8 

Ni-60 

Co-60m 

4 

13 

4 

Ni-58 

Co-58 

2 

2 

2 

Ni-58 

Fe-55 

7 

12 

6 

Fe-54 

Cr-51 

13 

20 

13 

Cr-53 

V-52 

37 

0 

0 

Cr-53 

V-53 

15 

25 

18 

Fe-54 

Mn-54 

6 

8 

7 

Cr-50 

Cr-51 

22 

3 

30 

Fe-57 

Mn-56 

36 

0 

0 

Fe-57 

Mn-57 

18 

26 

33 

Ni-62 

Fe-59 

19 

35 

71 

Ni-62 

Co-62m 

30 

50 

56 

Ni-62 

Co-61 

39 

79 

66 

Ni-62 

Co-62 

28 

46 

47 

Ni-60 

Co-60 

16 

28 

14 

Cr-54 

Cr-55 

24 

22 

67 

Cr-54 

Ti-51 

27 

49 

26 

Cr-54 

V-54 

25 

52 

50 

Fe-54 

Fe-55 

26 

14 

48 
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Figure 2: Comparison of the Alloy test case results, showing the first ten interactions 
according to the Pathways Based Metric Spbm , ranked accordingly. 
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3.2.2 Alloy+c 

See Table |3} Table [5j Table [6j Figure [3j Figure [4] and Figure [5} 

Table [4] suggests that once the Pearson correlation becomes below 0.1 it becomes inac¬ 
curate. Figure [3] shows that the smaller Pearson coefficients vary erratically with sampling, 
from which it is inferred that round-off effects have become important. 

As indicated in Table [l] this case involves both an irradiation phase and a cooling period. 
Care is required in comparing the FDM approach in this instance, for the method uses 
only the matrix for the cooling phase, whereas the other analyses are of the entire history. 
Although there is still reasonably good correlation between Pfds and Pprs , it is not as good 
in the other test cases. 



Figure 3: Comparison of the Alloy+c test case results, showing the first ten odd-numbered 
interactions according to the scaled Pearson technique value Sprs, for a Monte-Carlo sample 
size of N x = 640 per reaction, as N x is increased. 
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Table 4: Alloy+c case. Values of absolute Pearson correlation coefficient |r| as Monte Carlo 
sample size increases with N x . 




Absolute Pearson 

Parent 

Child 

40 

160 

640 

Ni-58 

Fe-55 

0.81342 

0.79097 

0.79709 

Fe-54 

Mn-54 

0.48716 

0.48798 

0.48761 

Ni-58 

Co-57 

0.24152 

0.23469 

0.23442 

Mn-55 

Mn-54 

0.17194 

0.16534 

0.16563 

Ni-60 

Co-60m 

0.16462 

0.15008 

0.13389 

Ni-58 

Co-58 

0.13754 

0.12544 

0.11864 

Fe-56 

Fe-55 

9.410 x 10" 2 

7.202 x 10" 2 

9.237 x 10" 2 

Ni-60 

Co-60 

0.10390 

6.218 x 10" 2 

5.669 x 10" 2 

Ni-58 

Co-58m 

6.996 x 10~ 2 

5.098 x 10" 2 

4.101 x 10~ 2 

Ti-46 

Sc-46 m 

3.614 x 10” 2 

1.570 x 10" 3 

1.575 x 10" 2 

Ti-47 

Sc-46 

2.538 x 10~ 2 

2.806 x 10~ 3 

1.209 x 10“ 2 

V-49 

Sc-46 

4.678 x 10" 3 

1.198 x lO’ 2 

1.144 x 10“ 2 

Co-60m 

Co-60 

1.923 x 10“ 2 

9.042 x 10~ 3 

8.303 x 10“ 3 

Cr-50 

V-50 

9.279 x 10~ 3 

1.598 x 10- 2 

8.073 x 10“ 3 

Ni-60 

Co-59 

1.648 x 10~ 2 

1.167 x 10~ 2 

7.639 x 10“ 3 

Ni-57 

Co-57 

1.665 x 10~ 2 

1.499 x lO’ 2 

7.389 x 10~ 3 

Co-57 

Co-58 

2.483 x 10~ 3 

1.572 x 10" 2 

6.780 x 10" 3 

Fe-54 

Cr-51 

2.407 x 10~ 2 

2.252 x 10" 2 

6.773 x 10" 3 

Co-57 

Co-58m 

3.906 x 10~ 2 

9.176 x 10~ 3 

6.739 x 10" 3 

Ti-47 

Ti-46 

2.828 x 10~ 2 

1.220 x lO’ 2 

5.516 x 10" 3 

Co-59 

Fe-59 

8.709 x lO’ 3 

1.451 x lO’ 2 

5.499 x 10“ 3 

Fe-58 

Fe-59 

2.626 x 10~ 2 

1.466 x 10- 2 

5.289 x 10“ 3 

Fe-54 

Fe-55 

2.698 x 10~ 2 

7.802 x 10~ 3 

5.261 x 10“ 3 

Co-58 

Fe-58 

2.502 x 10~ 3 

1.069 x 10" 3 

5.199 x 10“ 3 

Ti-47 

Sc-46 m 

4.156 x 10~ 3 

1.867 x lO’ 2 

5.039 x 10“ 3 

Cr-50 

Cr-51 

1.015 x 10“ 2 

2.073 x 10" 3 

4.815 x 10" 3 

Co-59 

Co-60 

1.789 x 10“ 2 

1.673 x lO” 2 

4.743 x 10" 3 

Ni-58 

Ni-59 

2.777 x 10" 3 

1.273 x lO” 2 

4.581 x 10" 3 

Mn-55 

Mn-56 

1.473 x 10~ 2 

3.741 x 10" 3 

4.456 x 10" 3 

Co-58 

Co-59 

5.478 x 10" 3 

6.033 x 10~ 4 

4.220 x 10" 3 
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Tabic 5: Alloy+c case. Rankings for different methods. 




Sensitivity 

Parent 

Child 

FDS 

PBM 

PRS 

Cr-52 

Cr-51 

11 

39 

49 

Mn-55 

Mn-54 

6 

1 

4 

Fe-56 

Fe-55 

9 

6 

7 

Ni-58 

Co-57 

5 

3 

3 

Ni-58 

Ni-57 

16 

20 

42 

Ni-58 

Co-58 

1 

7 

6 

Ni-58 

Co-58m 

2 

9 

9 

Ni-58 

Fe-55 

4 

4 

1 

Fe-54 

Cr-51 

8 

40 

17 

Fe-54 

Mn-54 

3 

2 

2 

Cr-50 

Cr-51 

13 

26 

26 

Ni-62 

Fe-59 

12 

37 

44 

Ni-60 

Co-60 

10 

14 

8 

Ni-60 

Co-60m 

7 

11 

5 

Fe-54 

Fe-55 

15 

5 

23 

Cr-50 

V-49 

14 

16 

35 

Fe-58 

Fe-59 

18 

27 

22 

Fe-54 

Mn-53 

17 

18 

30 

Ni-60 

Co-59 

22 

28 

15 

Ni-59 

Co-58 

21 

0 

0 

Ni-58 

Fe-54 

31 

0 

0 

Ni-59 

Co-58m 

23 

0 

0 

Co-59 

Fe-59 

19 

42 

20 

Co-59 

Co-58 

26 

0 

0 

Fe-56 

Mn-55 

30 

0 

0 

Co-59 

Co-58m 

25 

0 

0 

Ni-59 

Fe-55 

32 

0 

0 

Ni-62 

Ni-63 

27 

12 

33 

Fe-55 

Mn-54 

28 

0 

0 

Co-58 

Co-57 

24 

30 

34 

Co-57 

Co-56 

29 

44 

37 

Co-58 

Mn-54 

33 

0 

0 

Co-58 

Co-58m 

20 

62 

0 
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Figure 4: Comparison of the Alloy+c test case results, showing the first ten interactions 
by magnitude of Frechet derivative. The labels are ordered according to Frechet derivative 
size, so that the top interaction is the most sensitive. 



Figure 5: Comparison of the Alloy+c test case results, showing the first ten interactions 
according to the Pathways Based Metric Spbm for which comparison is possible. 
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Table 6: Alloy+c case. Rankings for different methods. 




Sensitivity 

Parent 

Child 

FDS 

PBM 

PRS 

Mn-55 

Mn-54 

6 

1 

4 

Fe-54 

Mn-54 

3 

2 

2 

Ni-58 

Co-57 

5 

3 

3 

Ni-58 

Fe-55 

4 

4 

1 

Fe-54 

Fe-55 

15 

5 

23 

Fe-56 

Fe-55 

9 

6 

7 

Ni-58 

Co-58 

1 

7 

6 

Co-60m 

Co-60 

0 

8 

13 

Ni-58 

Co-58m 

2 

9 

9 

Co-58m 

Co-58 

0 

10 

32 

Ni-60 

Co-60m 

7 

11 

5 

Ni-62 

Ni-63 

27 

12 

33 

Co-58 

Co-59 

0 

13 

31 

Ni-60 

Co-60 

10 

14 

8 

Co-59 

Co-60m 

0 

15 

43 

Cr-50 

V-49 

14 

16 

35 

Co-59 

Co-60 

0 

17 

27 

Fe-54 

Mn-53 

17 

18 

30 

Mn-53 

Mn-54 

0 

19 

46 

Ni-57 

Co-57 

0 

21 

16 

Ni-58 

Ni-57 

16 

20 

42 

Mn-55 

Mn-56 

0 

23 

29 

Mn-56 

Fe-56 

0 

22 

0 

Co-58m 

Co-59 

0 

24 

36 

Co-57 

Co-58 

0 

25 

18 

Cr-50 

Cr-51 

13 

26 

26 

Fe-58 

Fe-59 

18 

27 

22 

Ni-60 

Co-59 

22 

28 

15 

Ni-64 

Ni-63 

0 

29 

39 

Co-58 

Co-57 

24 

30 

34 
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3.2.3 Fe 


See Table [7] and Figure [6} 



Figure 6: Comparison of the Fe test case results, showing the first ten interactions according 
to the Pathways Based Metric Spbm, ranked accordingly. 
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Table 7: Fe case. Rankings for different methods. 




Sensitivity 

Parent 

Child 

FDS 

PBM 

PRS 

Fe-56 

Mn-56 

1 

1 

1 

Fe-56 

Fe-55 

4 

4 

3 

Fe-54 

Cr-51 

3 

5 

4 

Fe-54 

Mn-54 

2 

2 

2 

Fe-57 

Mn-56 

13 

0 

0 

Fe-57 

Mn-57 

5 

7 

6 

Fe-54 

Fe-55 

6 

3 

27 

Fe-56 

Mn-55 

7 

27 

18 

Fe-58 

Fe-59 

9 

6 

25 

Fe-58 

Cr-55 

10 

19 

8 

Fe-58 

Mn-58m 

11 

21 

19 

Fe-58 

Mn-58 

12 

22 

22 

Fe-54 

Mn-53 

8 

9 

7 

Fe-56 

Fe-57 

16 

10 

11 

Fe-55 

Mn-54 

14 

0 

0 

Cr-54 

Cr-55 

20 

15 

12 

Cr-54 

V-54 

21 

26 

5 

Mn-55 

Cr-55 

18 

24 

26 

Mn-55 

V-52 

22 

25 

16 

Mn-55 

Mn-56 

15 

0 

0 

Fe-57 

Fe-58 

25 

13 

24 

Mn-55 

Mn-54 

24 

0 

0 

Fe-57 

Cr-54 

27 

28 

14 

V-51 

V-52 

32 

16 

23 

Mn-53 

Mn-54 

19 

8 

10 

Fe-57 

Fe-56 

23 

0 

0 

Co-59 

Co-60m 

30 

12 

21 

Co-59 

Mn-56 

31 

0 

0 

Co-59 

Fe-59 

28 

0 

0 

Fe-55 

Mn-55 

17 

23 

17 
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3.2.4 LiMix 


See Table [8] and Figure [7} The comparison between the various metrics in Figure [7] does 
not at first appear to be as successful as in other cases. However the dominant interaction 
from the PBM involves tritium for which uncertainty data are not accessible in the database, 
hence the FDS and PRS cannot assign it a ranking and it is omitted from the plot. Moreover 
all FDS rankings over 21 similarly correspond to zero uncertainty and allowing for this, the 
comparison is as good as any reported herein. 



Figure 7: Comparison of the LiMix test case results, showing the first ten interactions 
according to the Pathways Based Metric Spbm for which comparison is possible. 
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3.2.5 WMix 


See Table |9j , Table [TOj Table |TT| Figure [8j Figure [9] and Figure p~0| Table [9] suggests that 
once the Pearson correlation becomes below 0.1 it becomes inaccurate. Figure [8] shows that 
the lower rankings in terms of sensitivity vary erratically with sampling for similar reasons 
to do with round-off effects. 



Figure 8: Comparison of the WMix test case results, showing the first ten odd-numbered 
interactions according to the scaled Pearson technique value Sprs, for a Monte-Carlo sample 
size of N x = 640 per reaction, as N x is increased. 
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Table 9: WMix case. Values of absolute Pearson correlation coefficient |r| as Monte Carlo 
sample size increases with N x . 




Absolute Pearson 

Parent 

Child 

40 

160 

640 

Rel87 

Rel88 

0.93702 

0.94033 

0.94091 

Irl93 

Irl94 

0.24971 

0.23407 

0.24135 

Irl91 

Irl92 

0.13306 

0.13520 

0.13542 

Irl93 

Irl93m 

0.11155 

0.11969 

0.11550 

Re 185 

Rel86 

0.10156 

8.413 x 10~ 2 

7.955 x 10~ 2 

W-184 

W-185 

6.398 x 10“ 2 

9.215 x 10" 2 

7.856 x 10" 2 

W-186 

W-187 

5.029 x 10“ 2 

6.135 x 10- 2 

6.423 x 10" 2 

Re 187 

Rel88m 

2.181 x 10~ 2 

5.108 x 10" 2 

5.264 x 10" 2 

Ptl92 

Ptl91 

4.270 x 10“ 2 

4.238 x 10" 2 

3.215 x 10" 2 

Irl93m 

Irl 93 

2.090 x 10" 3 

9.212 x 10" 3 

2.514 x 10" 2 

Irl91 

Irl92m 

3.290 x 10~ 2 

1.785 x 10~ 2 

2.510 x 10“ 2 

W-186 

W-185m 

4.033 x 10“ 2 

4.676 x 10- 2 

2.323 x 10“ 2 

Re 187 

Rel86 

3.941 x 10~ 2 

2.248 x 10~ 2 

1.720 x 10“ 2 

Bi209 

Bi210 

6.132 x 10“ 5 

1.141 x 10- 2 

1.679 x 10“ 2 

Irl92 

Irl 93m 

1.539 x 10~ 2 

1.425 x 10- 2 

1.637 x 10“ 2 

Irl92 

Irl 93 

1.845 x 10~ 2 

1.529 x 10- 2 

1.352 x 10" 2 

W-182 

W-181 

1.497 x 10~ 2 

1.685 x 10" 2 

1.274 x 10" 2 

Irl91 

Irl 90 

2.775 x 10~ 2 

8.743 x 10~ 3 

1.263 x 10" 2 

Pb208 

Pb207m 

3.345 x 10~ 2 

9.360 x 10~ 3 

1.217 x 10" 2 

Irl91 

Irl 91m 

2.715 x 10~ 2 

1.390 x 10" 2 

1.167 x 10" 2 

Ptl94 

Ptl93m 

1.482 x 10~ 2 

1.318 x 10- 2 

1.020 x 10“ 2 

W-186 

W-185 

3.379 x 10~ 2 

1.427 x 10- 2 

9.663 x 10“ 3 

W-183 

W-183m 

2.160 x 10~ 3 

1.927 x 10- 2 

9.005 x 10“ 3 

W-183 

W-184 

5.385 x 10- 3 

1.532 x 10~ 3 

8.812 x 10“ 3 

Irl93 

Irl 92 

2.121 x 10~ 2 

8.499 x 10~ 3 

8.497 x 10“ 3 

Re 185 

Rel84m 

1.874 x 10“ 2 

7.629 x 10" 3 

8.467 x 10" 3 

Irl94 

Irl 95m 

4.688 x 10" 3 

5.884 x 10" 3 

8.447 x 10" 3 

Irl93 

Osl93 

3.587 x 10“ 4 

2.355 x 10" 3 

8.253 x 10" 3 

Re 188m 

Rel88 

1.010 x 10~ 2 

4.477 x 10" 3 

7.839 x 10" 3 

Bi210m 

Bi210 

1.836 x 10" 2 

3.068 x 10~ 3 

7.527 x 10" 3 
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Table 10: WMix case. Rankings for different methods. 




Sensitivity 

Parent 

Child 

FDS 

PBM 

PRS 

Bi209 

Bi210 

14 

38 

14 

Re 187 

W-185m 

29 

0 

0 

Re 187 

Rel88m 

6 

12 

8 

Irl93 

Osl91 

35 

0 

0 

Irl93 

Osl91m 

37 

0 

0 

Bi209 

Pb207m 

28 

0 

0 

Irl93 

Irl92m 

11 

0 

0 

Re 187 

W-185 

31 

0 

0 

W-184 

W-185m 

17 

50 

51 

Re 187 

W-187 

19 

0 

0 

Irl93 

Irl 93m 

3 

16 

4 

Re 187 

Tal83 

32 

0 

0 

W-186 

W-185m 

7 

27 

12 

Re 187 

Rel86 

8 

18 

13 

Re 187 

Rel88 

1 

4 

1 

Ptl94 

Osl91 

26 

0 

0 

Ptl94 

Osl91m 

30 

0 

0 

Pb208 

Pb207m 

12 

26 

19 

Re 185 

W-185m 

20 

0 

0 

Irl93 

Irl 92 

9 

20 

26 

W-184 

W-185 

4 

13 

6 

Irl93 

Irl 94 

2 

1 

2 

Ptl94 

Irl 93m 

34 

0 

0 

W-186 

W-185 

10 

21 

22 

W-184 

Tal83 

33 

0 

0 

W-184 

Tal82m 

38 

0 

0 

W-184 

Tal82 

36 

0 

0 

W-184 

W-183m 

15 

35 

55 

Ptl94 

Irl 94 

24 

0 

0 

Ptl94 

Ptl93m 

13 

24 

21 
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Figure 9: Comparison of the WMix test case results, showing the first ten interactions by 
magnitude of Frechet derivative. The labels are ordered according to Frechet derivative size, 
so that the top interaction is the most sensitive. 



Figure 10: Comparison of the WMix test case results, showing the first ten interactions 
according to the Pathways Based Metric Spbm for which comparison is possible. 
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Table 11: WMix case. Rankings for different methods. 




Sensitivity 

Parent 

Child 

FDS 

PBM 

PRS 

Irl93 

Irl94 

2 

1 

2 

Re 185 

Rel86 

0 

2 

5 

Irl91 

Irl92m 

0 

3 

11 

Re 187 

Rel88 

1 

4 

1 

Irl91 

Irl92 

0 

5 

3 

Irl92m 

Irl92 

0 

6 

71 

lrl92 

Irl93m 

0 

7 

15 

W-186 

W-187 

5 

8 

7 

Irl93m 

Irl 93 

0 

9 

10 

Irl92 

Irl 93 

0 

10 

16 

W-187 

Rel87 

0 

11 

0 

Re 187 

Rel88m 

6 

12 

8 

W-184 

W-185 

4 

13 

6 

Re 186 

W-186 

0 

14 

37 

Re 188m 

Rel88 

0 

15 

29 

Irl93 

Irl 93m 

3 

16 

4 

W-185 

Rel85 

0 

17 

0 

Re 187 

Rel86 

8 

18 

13 

Irl94 

Irl 95 

0 

19 

62 

Irl93 

Irl 92 

9 

20 

26 

W-186 

W-185 

10 

21 

22 

W-182 

W-181 

0 

22 

17 

Pb207 

Pb207m 

0 

23 

54 

Ptl94 

Ptl93m 

13 

24 

21 

Irl94 

Ptl94 

0 

25 

0 

Pb208 

Pb207m 

12 

26 

19 

W-186 

W-185m 

7 

27 

12 

Irl91 

Irl 91m 

0 

28 

20 

Re 186 

Osl86 

0 

29 

0 

Osl86 

Osl85 

0 

30 

38 
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Table 12: Y203 case. First five radio-nuclides in final inventory, ordered by activity, plus 


Sr-89. H 

alf-lives from the EASY-1I database. 

Order 

Nuclide 

Activity 

Atoms 

Half-life 



Bq 

Percent 


Tk 

1 

Y-89m 

2.5874 x 10 7 

89.28% 

5.847 x 10 8 

15.663s 

2 

N-16 

2.9911 x 10 6 

10.32% 

3.077 x 10 7 

7.13s 

3 

Rb-86m 

1.0193 x 10 5 

0.35% 

8.971 x 10 6 

1.017m 

4 

C-15 

1.0083 x 10 4 

0.04% 

3.562 x 10 4 

2.449s 

5 

Y-88 

1.191 x 10 3 

0.004% 

1.583 x 10 10 

107d 

9 

Sr-89 

62 

0.0002% 

3.891 x 10 s 

50.57d 


3.2.6 Y203 

This very simple case illustrates the ranking process in additional detail, giving examples of 
the nuclear data that is used in the calculations. Following irradiation, the activity at the 
end of this test case is dominated (99.6 %) by two nuclides N-16 and Y-89m, see Table [l2| for 
details. The pathways analysis performed routinely by FlSPACT-II shows that all pathways 
leading the 6 nuclides listed, consist of just one reaction. 


PRS As a result of the dominance by two reactions, all the detailed rankings by Pearson 
apart from the first two are suspect, see Table [14| and Figure pT] in support of this contention. 


Moreover, no (zero) uncertainty estimate is provided for the 0-161 N-16 reaction, see Table 13 


the source of the N-16 in the final inventory. Thus the Pearson ranking is maximal, and in 
the other sections 0-16|N-16 would have to be omitted from the comparison plots. Other 
important reactions are identified by inspection of comparison tables and plots such as 
Table 15 and Figure [Irrespectively. On this basis, the reaction producing Sr-89 is identified 


as potentially important, hence its inclusion in Table 12 
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Table 13: Y203 case. Reactions identified as important by pathways analysis. Note that 
apart from the reactions involving an excited parent species (labelled ‘m’) each reaction 
is equivalent to a pathway. For the notation describing each reaction see the Fispact-II 
manual._ 


Pathway 

Reaction 

Cross-section 

barn 

€ij Relative 
Uncertainty % 

0-16(n,p)N-16 

0.03357 

0.0 

0-17(n,p)N-17 

0.01020 

27.0 

0-18(n,a)C-15 

0.05645 

28.0 

0-18(n,d)N-17 

8.9740 x 10" 6 

23.0 

Rb-86m(b)Rb-86 

2.2675 

4.39 

Y-89(n,a)Rb-86 

0.0036110 

70.0 

Y-89(n,a)Rb-86m 

0.001820 

70.0 

Y-89(n,p)Sr-89 

0.02179 

19.0 

Y-89(n,2n)Y-88 

0.88630 

5.7 

Y-89(n,n)Y-89m 

0.4347 

8.5 

Y-89(n,g)Y-90 

0.0022040 

9.2 

Y-89(n,g)Y-90m 

2.1410 x 10” 4 

9.2 

Y-90m(n,g)Y-90 

2.4633 

4.51 


Table 14: Y203 case. Values of absolute Pearson correlation coefficient |r| as Monte Carlo 
sample size increases with N x . 




Absolute Pearson 

Parent 

Child 

160 

640 

2560 

Y-89 

Y-89m 

0.99945 

0.99945 

0.99948 

Y-89 

Rb-86m 

5.664 x 10“ 2 

4.483 x 10“ 2 

2.851 x 10“ 2 

0-18 

C-15 

1.768 x 10" 3 

1.504 x 10“ 4 

8.698 x 10" 3 

Y-90m 

Y-90 

1.041 x 10" 2 

1.250 x 10" 2 

7.383 x 10" 3 

Y-89 

Rb-86 

1.326 x 10" 3 

1.010 x 10" 2 

3.955 x 10" 3 

0-18 

N-17 

2.644 x 10" 2 

8.750 x 10" 3 

3.389 x 10" 3 

Y-89 

Y-90 

3.402 x 10" 3 

6.061 x 10" 3 

2.537 x 10" 3 

Y-89 

Y-90m 

1.537 x 10" 2 

2.666 x 10" 3 

2.261 x 10" 3 

Rb-86m 

Rb-86 

6.646 x 10 -3 

1.651 x 10“ 2 

2.238 x 10" 3 

0-17 

N-17 

1.821 x 10“ 2 

2.273 x 10“ 2 

1.766 x 10“ 3 

Y-89 

Y-88 

2.699 x 10“ 2 

1.402 x 10“ 2 

1.426 x 10“ 3 

Y-89 

Sr-89 

1.358 x 10“ 2 

1.085 x 10“ 2 

4.010 x 10 -4 

0-16 

N-16 

2.585 x 10~ 16 

5.905 x 10~ 16 

1.552 x 10“ 15 
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CVJ 


Figure 11: Comparison of the Y203 test case results, showing the first ten interactions 
according to the scaled Pearson technique value Spps, f° r a Monte-Carlo sample size of N x = 
2 560 per reaction, as N x is increased. 









FDS The Frechet derivative method works with a matrix that has coefficients of all the 
possible reactions which involve the species by the standard Fispact-II pathways analysis. 
There seems to be little point in listing them all, obviously Table 13 is indicative. As in 
the case of Pearson, the absence of an uncertainty estimate for the 0-16|N-16 reaction leads 
to a maximal ranking for the scaled Frechet derivative, and in the other sections 0-16|N-16 
would have to be omitted from the comparison plots. 



Figure 12: Comparison of the Y203 test case results, showing the first ten interactions by 
magnitude of Frechet derivative. The labels are ordered according to Frechet derivative size, 
so that the top interaction is the most sensitive. 
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Table 16: Y203 case. Rankings for different methods. 




Sensitivity 

Parent 

Child 

FDS 

PBM 

PRS 

Y-89 

Y-89m 

2 

1 

1 

0-16 

N-16 

30 

2 

13 

Y-89 

Rb-86m 

4 

3 

2 

0-18 

C-15 

7 

4 

3 

Y-89 

Y-88 

1 

5 

6 

0-17 

N-17 

10 

6 

7 

Y-89 

Y-90m 

8 

7 

11 

Y-89 

Y-90 

6 

8 

10 

Y-89 

Sr-89 

3 

9 

9 

Y-89 

Rb-86 

5 

10 

5 

0-18 

N-17 

12 

11 

12 

Rb-86m 

Rb-86 

0 

12 

8 

Y-90m 

Y-90 

0 

13 

4 


PBM Since all pathways leading to the most important nuclides in terms of activity, 
consist of just one reaction, it follows that these reactions are simply ranked in order of 
contribution of the child species to the final activity. Even though it has no associated 
uncertainty, the 0-16|N-16 reaction can be ranked by PBM. Taking into account the fact 
that in the other sections 0-16|N-16 would have to be omitted from the comparison plots, 
Figure [13] indicates that there is surprisingly good agreement between the methods, even for 
reactions that contribute little to the final activity. 
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Figure 13: Comparison of the Y203 test case results, showing the first ten interactions 
according to the Pathways Based Metric Spbm- 


32 









4 Conclusions 


The sensitivity of the total activity of an inventory to uncertainties in the nuclear data 
for neutron-induced reactions has been studied. Six different test cases covering nearly the 
whole range of atomic masses were considered using three different ranking techniques. It 
is expected that similar results would be obtained for other inventory properties and other 
projectile particle species. 

The principal result is that a simple pathways-based metric (PBM) gives a sensitivity 
ranking of interactions which is comparable to ranking based on more conventional mea¬ 
sures obtained either by the direct method or in terms of Pearson correlation coefficients. 
Moreover, the PBM is superior in that it 

1. is quick to calculate once the principal pathways have been identified 

2. does not suffer from numerical difficulties such as underflow (Frechet direct) or round¬ 
off (Pearson) in its evaluation 

3. may be generalised to the case of multiple irradiation periods just like the pathways- 
reduced approach itself, whereas the other two techniques require further investigation. 

4. does not require error estimates for every interaction coefficient like Pearson. 

An additional noteworthy feature is that the PBM, which is a global measure of un¬ 
certainty, is comparable with more local measures, provided these others are scaled by the 
uncertainty in the reaction cross-section. This scaling is to be expected since the uncer¬ 
tainty estimates computed by Fispact-II |T, §A.13] involve a multiplication by a measure 
of cross-section uncertainty (r.m.s. is used to combine reaction coefficients rather than the 
simple percentages). However, the product also involves the number of child nuclides in the 
inventory which is a significantly different measure from the point sensitivity measures. 

The value of studying a wide range of test cases is that it demonstrates the general 
applicability of the above conclusions. In conjunction with modifications to FiSPACT-II for 
more efficient pathways-based analysis in the presence of multiple irradiations, the PBM 
should be extended to account for loops in the pathways and ultimately integrated into a 
production version of the Fispact-II package. 
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Appendix: Frechet Derivatives 

As explained in Section [TJ the Bateman equation Eq. ([Tj) 

dX 

— = ax, x(o) = x Q , te[o,t f ], 

where X G M n is a vector of nuclide numbers and A G M nxn is a matrix of nuclear interaction 
coefficients, controls the evolution of the nuclear activation over time. In this appendix, we 
focus on the case where A is constant in time. 

We are interested in the sensitivity of the total activity Eq. (|2]) 

n 

Q = ^ 4 x k {tf) 

k=1 

to the elements in A, which is determined by the n 2 numbers dQ/dAij. To determine these 
quantities we use the matrix exponential and its Frechet derivative. The matrix exponential 
of A G M nxn is defined by 



The Frechet derivative of the exponential at A in the direction E G M. nxn is denoted by 
L exp (A, E) G M nxn and satisfies 

e A+E =e A + Lexp(A! E)+0 (|| E ||). 

For further details of Frechet derivatives see jUl Chap. 3]. 

The solution to the Bateman equation is given by 

X(t) = e M X 0 


and so 

Q = f T X(t f ) = f T e Mf X 0 , f = [h . .. l n ] T . 

Let E ij be the n x n matrix with a 1 in the (i,j) entry and zeros elsewhere. Now, 

dQ _ Q(Aij + h) — Q(Aij) 
dAij s—S 

fT (g(A+E VJ 8)t f _ e A tf\ Xo 

= lim------- 

< 5->0 0 

= lim f T (L exp (Af / ,E u V) + O (h))X 0 

< 5—>-0 S 

= t ff J -^exp (At/, Ey)X 0 , 

where we have used the fact that L exp is linear in its second argument. 


34 






To determine the k largest of these derivatives we can simply compute them all and sort 
them. For this we can use the relationship m eq. (3.16)] 


( 

'tA 

LU 


e tA L exp (fA, E i:j ) 

l 

i 

O 

tA 


0 e tA 


which yields the formula 


( 

'tA 

Ep- 


' 0 ' 


-^exp 

l 

0 

tA 

) 

A 0 _ 


< 

X 

_ i 


(15) 


(16) 


Hence one method to compute L exp (Atf, Eij)X 0 is to apply the method from [HI] to compute 
the product on the left-hand side and then read off the first n components. 

However, it is not necessary to carry out n 2 Frechet derivative evaluations. One suffices, 
as we now explain. We need some notation. The Kronecker product of two matrices is now 
text-book for numerical linear algebra, see e.g. ref HU. It is defined for matrices B and C 
(of any dimension) as the block matrix B ® C = (b tJ C). The vec operator stacks the columns 
of a matrix one of top of each other from first to last, producing a long vector. We need 
the property that vec(L exp (A, E)) = K(A) vec(E), for some n 2 x n 2 matrix K(A) that satisfies 
K(A) 7 = K(A T ). Using the fact that the vec of a scalar is itself and the formula 


vec(AXB) = (B 7 ® A)vec(X), 


we have 


dQ 

dA i:j 


tff T L exp (At,f, E ij)X, 


o 


vec (tff T L exp (At f , Eij)X 0 ) 
tf(X o ® f) T vec (L exp (Af / , Ey)) 
^/5 |T K(Af/) vec(Ejj), 


where g = Xq ® /. Now, since vec(Ey) is a unit vector, we simply require the k largest 
elements in modulus of g 1 K(A tf), which are the largest k elements in magnitude of K(Atf) T g. 
We have K(A T tf)g = vec(L exp (A 7 tf, E)), where vec(E) = g = X 0 ® f and hence E = fx p . 
This means that a single Frechet derivative evaluation is sufficient, and it can be done using 
the relationship (15) above with an algorithm to compute the matrix exponential such as 
that in [ IS] . 

The computation of the exponential requires numerous matrix products, which can oc¬ 
casionally cause numerical over- or under-flow due to the large range of magnitudes in the 
coefficients arising in nuclear activation problems. This may necessitate the use of quadruple 
precision arithmetic on certain problems. In fact, quadruple precision was used to check the 
accuracy of all the Frechet derivatives calculated in the course of the current work. 
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